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Abstract. Schwarzschild's solution to the Einstein Field Equations was one of the first and 
most important solutions that lead to the understanding and important experimental tests 
of Einstein's theory of General Relativity. However, Schwarzschild's solution is essentially 
based on an ideal theory of gravitation, where all inhomogeneities are ignored. There- 
fore, any generalization of the Schwarzschild solution should take into account the effects 
of small perturbations that may be present in the gravitational field. The theory of Macro- 
scopic Gravity characterizes the effects of the inhomogeneities through a non-perturbative 
and covariant averaging procedure. With similar assumptions on the geometry and matter 
content, a solution to the averaged field equations as dictated by Macroscopic Gravity are 
derived. The resulting solution provides a possible explanation for the flattening of galactic 
rotation curves, illustrating that Dark Matter is not real but may only be the result of 
averaging inhomogeneities in a spherically symmetric background. 



1. Introduction 

The most common assumption in the study of modern cosmology and gravity is that the 
gravitational field is described by Einstein's Theory of General Relativity which is assumed 
to be valid on all length scales. Cosmologies and isolated gravitational models based on 
this assumption have been very successful in explaining a wide range of observational phe- 
nomenon. For example, on length scales of the solar system, the bending of light curves, 
the precession of perihelia, the gravitational red-shift, and the radar echo delay all provide 
significant evidence that Einstein's Theory of General Relativity is adequate to describe the 
gravitational dynamics in our own solar system pp. Furthermore, the observed expansion of 
the universe, the observed abundance of light elements created in the first few moments after 
the Big Bang, and the observed relic Cosmic Microwave Background radiation all attest to 
the validity of Einstein's theory of General Relativity on much larger scales. 

However, everything is not as rosy as it might seem. Contrary to Einstein's theory, the 
universe is accelerating [2j. This acceleration cannot be explained in General Relativity 
without the introduction a mysterious substance labeled "Dark Energy" , a key ingredient in 
the popular and widely accepted ACDM model of the universe. The existence and properties 
of Dark Energy are still a matter of debate, but a different solution may be more attractive, 
that is, perhaps an alternative theory of gravity ought to be used when modeling large 
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scale gravitational systems such as the universe. The question then becomes, can General 
Relativity provide an accurate description of the gravitational field on all length scales? 

Indeed, even on galactic scales, both Newtonian gravity and Einstein's theory of General 
Relativity struggle to resolve the observed rotation curves of spiral galaxies[3j. There is 
currently no way to explain the flatness of the rotation curves as the distance from the galactic 
center increases without introducing a substance that has been named "Dark Matter" , the 
second ingredient in the popular ACDM model. Dark Matter is a novel and unknown 
material that does not emit any light, and therefore is difficult to observe directly. However, 
even the premise of dark matter has its problems; it is sometimes modeled as a cloud of 
pressureless dust, and therefore, the mass distribution of the dark matter would behave as 
r~ 3 , whereas a solution to the flatness of the rotation curves would really require a mass 
distribution having a profile of r~ 2 . Again, perhaps an alternative and more natural solution 
exists which explains the dynamics of gravitational systems on intermediate scales, i.e., on 
the scale of galaxies. 

Both the universe and a galaxy can be loosely described as systems of point-like or discrete 
particles. For the universe, the constituent particles would be galaxies or galaxy clusters 
while for a galaxy the constituent particles would be the stellar systems within. Because we 
are dealing with inhomogeneous distributions of point-like particles in each case, the problem 
of determining the gravitational dynamics (inside the universe or inside a galaxy) is one that 
requires either a relativistic N-body simulation or an averaging procedure that averages both 
sides of the Einstein Field equations on some appropriate scale. 

Given that General Relativity has been tested to a high degree of accuracy on scales 
of the solar system P, it is therefore reasonable to assume that any averaging procedure 
used to model large scale gravitational effects ought to be based on Einstein's theory of 
General Relativity. There have been several approaches that purport to average gravitational 
fields based on Einstein's Theory of General Relativity, (see [U El El El E] and references 
within for a complete discussion). However, almost all of these approaches are either non- 
covariant, or perturbative in nature and therefore their general applicability is questionable. 
The theory of Macroscopic Gravity as proposed and developed by R. Zalaletdinov [5j E] 
is covariant and non-perturbative by construction and therefore does not suffer the same 
weaknesses as the other approaches. The Zalaletdinov approach employs a covariant space- 
time averaging procedure for microscopic tensor fields that produces an averaged tensor 
field. When this averaging procedure is applied to the field equations of General Relativity, 
one obtains a generalization of the Einstein Field Equations of General Relativity for the 
averaged spacetime that includes a new tensor field. This new tensor field, constructed from 
a connection correlation tensor Z (which is further dependent upon an affine deformation 
tensor A) , satisfies an independent set of differential and algebraic equations [51 E] • 

Unfortunately for a variety of reasons including but not limited to the complexity of the 
equations describing Macroscopic Gravity, the potential of Macroscopic Gravity has not 
been fully explored; with the exception of [HI [101 [11] who have made preliminary attempts to 
understand the theory. Here, in this paper we will attempt to illustrate some of the effects of 
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the gravitational stress energy tensor due to the connection correlation field in a class of static 
spherically symmetric macroscopic spacetimes. The proposed gravitational models can be 
used to model galaxies in general, and galactic bulges and halos in particular. Reasonable and 
physically plausible assumptions on the macroscopic geometry and the connection correlation 
tensor will be made and the complete set of macroscopic gravity equations will be analyzed. 
The resulting solutions will not be solutions to the Einstein Field equations, but will be 
solutions to the Averaged Field Equations of Macroscopic Gravity. We will show that the 
effective energy density due to the connection correlation field has a mass distribution profile 
of r~ 2 , thereby providing an alternative explanation of the flattened galactic rotation curves. 

We shall use the same notation as [B] for all quantities. For quick reference, G a p is the 
metric of the averaged or macroscopic spacetime having a Levi-Cevita connection T p 7 = 
(J- a p.y), where T a ^ are the bilocal extensions of the Levi-Cevita connection of the micro- 
scopic spacetime. The Riemann curvature tensor corresponding to the macroscopic connec- 
tion T py is denoted as M a ^ with a Ricci tensor M Q/ g = aAt( g. Covariant differentiation 
with respect to connection T° ' p 1 is denoted with 1 1 and any index that is underlined is not in 
included in the antisymmetrization. The connection correlation tensor Z, is defined as the 
difference between the average of the product and the product of the averages of the bilocal 
extensions of the microscopic connection, via 

Please see papers [5] E] for a complete description of all equations and quantities involved. 
Given a coordinate system, and macroscopic metric G a p, many tensorial objects are greatly 
simplified by choosing an orthonormal vector bases or tetrad [T2], {e( Q )} (parentheses sur- 
rounding indices are used to denote tetrad indices) adapted to the metric such that the set 
of scalars representing the metric tensor 

are constant and take on values V(a){p) — diagfl, 1,1,-1]. Units are chosen so that c — G—1 
with k = 8n. 

2. The Macroscopic Gravity Equations 

2.1. Assumptions. Due to the difficulty in dealing with the extreme size of the tensorial 
objects involved in Macroscopic Gravity, we are encouraged to make reasonable geometric as- 
sumptions about all objects involved in order to make some progress into the understanding 
of the structure of the Macroscopic Gravity equations. One must recall that only when geo- 
metric assumptions were made about the spacetime (i.e., spherically symmetric and static) 
in General Relativity did one make progress in obtaining solutions of the Einstein Field 
Equations and insight into their physical meaning. 

Assumption 1: The Averaging Procedure and Length Scale. As outlined above we shall as- 
sume that the appropriate averaging procedure to use on large scales is the non-perturbative 
and covariant procedure as developed by Zalaletdinov [5] . In this paper, the averaging length 
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scale is assumed to be significantly larger than the size of our solar system, but much smaller 
than the size of typical galaxies (for example, on the order of 10-20 parsecs). 

Assumption 2: The Metric and Higher Order Connection Correlations. The metric correla- 
tions as outlined in [5] are assumed to be zero. In addition, we shall assume that the average 
of the inverse microscopic metric is equal to the inverse of macroscopic metric. This should 
be expected and is not considered unusual in any spacetime with a high degree of symmetry. 
Furthermore, all higher order connection correlations (except Z) of Macroscopic Gravity 
are assumed to be zero, this can be done in a self-consistent way as outlined in [5]. 

Assumption 3: Macroscopic Geometry. For the present moment we are interested in investi- 
gating the Macroscopic Gravity equations in static spherically symmetric spacetimes exterior 
and interior to some matter source. The assumption of spherical symmetry and staticity is 
consistent with the usual assumptions made about the geometry of the dark matter halo 
surrounding a galaxy. The macroscopic metric, G a p, can therefore be expressed without loss 
of generality, as 

(1) ds 2 = G aP dx a dx p = r 2 d6 2 + r 2 sin 2 6d<p 2 + e 2X{r) dr 2 - e 2u{r) dt 2 . 

The metric, equation (CQ), is invariant under a four-dimensional group of Killing vectors; 
the Killing vector, X = J^, and the three dimensional rotation group. This spacetime 
also allows one to define a timelike unit vector field, u a = [0, 0, 0, e~ u ^}, which has a non- 
zero acceleration, but has zero vorticity, zero shear, and zero expansion. The tetrad basis 
corresponding to metric ([1]) that will be used in future presentations is 

(2) e ( 3 = [i 0,0,0], = [0,^,0,0], 

e {3 ^ = [0,0,e- A W,0], e ( J = [0,0,0,e-^]. 

Assumption J^: Invariance of Macroscopic Gravity Objects. We shall assume that the con- 
nection correlation tensor Z is invariant under the same group of Killing vectors as the 
macroscopic metric. This assumption ensures that Z is compatible with the geometry of the 
macroscopic space. Note, Z need not be invariant under the same group of Killing vectors as 
the macroscopic spacetime a priori. In addition, we shall assume that the affine deformation 
tensor A has tetrad components that are functions of the coordinate r only. 

Assumption 5: Averaged Microscopic Matter. We shall also assume that the average of the 
microscopic stress energy tensor can be modeled as a perfect fluid, that is 

(3) (tj (micro) ) = p^u a u p + p^ {8% + u a u p ) 

where p mat and p mat are the energy density and pressure for the averaged matter and where 
u a can now also be interpreted as the average four-velocity of the fluid. We shall further 
assume that the energy density and the pressure of the averaged matter vanishes in the 
exterior of the mass producing the gravitational field. While in the interior of the mass 
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generating the gravitational field, the energy density of the averaged matter is assumed to 
be constant and non-zero, that is, 



(4) P 



(mat) 



p , r < R 

0, r> R 



where R is the size of the gravitational source. These assumptions on the averaged micro- 
scopic matter will yield the usual interior and exterior Schwarzschild solutions in the absence 
of connection correlations. 

Assumption 6: Electric part of the Connection Correlation Tensor. We shall assume that the 
electric part of the connection correlation tensor (see [13j E]) is zero, that is, ZE a (p ^ = 
Z a p i Ps v u v = 0. The primary motivation for setting ZE a j 3 1 $^ = is that the constraint 
equation fl5J) is identically satisfied. 

Z5 6 rva n | iy5 11 ryd a i rya 8 ry U 8 i 

f3[-y ktt^ 8e up] "I ^ /3[7 ua^ kit 8e] + ^ /% vc^ 8e ktt] + 

^CJ Z; $ e Zj K7T vcr] ~T Z' ^[7 ^cr ktt] T Z' /3[7 ktt^ 5e i/o"] — U. 

The assumptions 1 — 6 greatly simplify the extremely complex system of Macroscopic 
Gravity Equations, allowing us to analyze the effects of the gravitational stress-energy tensor 
(due to connection correlations) in the averaged gravitational field equations. 

2.2. The Connection Correlation Tensor. The connection correlation tensor Z, with 
coordinate components Z a j 3 tl 1 $ u , having index symmetries 

{ r}\ r/Oi 7 ryOi 7 77 Q 77 Q 

[V) Zj f3p, 8v — —A f3u Sfi — A 8v Pn — ~^ Sfi /3v 

can be uniquely represented by the set of scalars ^^ctHx/^MW) through the use of the 
tetrad ([2]). Assuming that Z is static and invariant under rotations places a significant 
number of constraints on the number of independent components. We find that after all 
symmetry conditions are imposed, only 150 independent components of Z remain, each of 
which is a function of the coordinate r only. With the assumption that ZE a p y 5 iJL = 0, we 
find that there are only 75 independent components of Z remaining. The algebraic cyclic 
identities, 

(7) Z a ^ Su] = 

yield an additional 50 linearly independent constraints. The differential cyclic constraint 
equations, 

(8) Z a p^ 1 Sy\\a] = 

can be sub-divided into two parts, an algebraic part (containing no derivatives), and a 
differential part. The algebraic part of the differential cyclic identities yields an additional 22 
algebraic constraints thereby leaving us with only 3 independent components of Z remaining. 
The differential part of the differential cyclic identities yields three differential equations for 
each of the three remaining independent variables. The remaining differential equations are 
easily solved with the result that the remaining independent tetrad components of Z behave 
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as ^^(<r)(x)^(K)(^)( r ) oc r 2 . Introducing three constants of integration hi,h 2 ,h^ the 32 
nontrivial components of Z are 

hi 





(2) 


Z (l) 


(D(i) (1) 


7(2) 


(2) 
(2)(2) V 




(1) 
(2)(2) V 


z® 


(3) 

(1)(1) 




(4) 
(1)(2) V 




(3) 
(1)(2) V 


Z (4) 


(3) 
(2)(2) V 



hi 

k 
hi 



where the remaining 24 nontrivial components are obtained via the symmetries of Z (see 
equation EJ). 

2.3. The Affine Deformation Tensor. The equations of Macroscopic Gravity also involve 
a set of equations that determines an affine deformation tensor A (see [5j |6] for additional 
details). The affine deformation tensor satisfies the algebraic equation 

(9) A t ( 3 [ p R a ^\] - A a e [ p R e p a x] = 0, 
and the differential equation 

( 10 ) A a p[ a \\ p ] - A a e [pA e p a ] = ~-Q a f3pa- 
where 

is a non-Riemannian curvature tensor and where 

Equation yields 32 constraints and therefore there are only 8 independent components of 
A remaining. The algebraic part equation (flUj) yields an additional 7 algebraic constraints 
and therefore there is only one independent component of A remaining. Solving the dif- 
ferential part of equation (TlOT) yields the last remaining component. The nontrivial tetrad 
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4(1) 

A (1)(1) 
AW, 



(2)(2) 
^ (2) (D( 2 ) 



r 

2^ 

r 

2^ 

r 



A( 4 ) (4)(4) = Ae^ (r) 
where Ao is a constant of integration. We observe that the sign of hi is now restricted, 
(11) h!>0. 



2.4. The Gravitational Stress-Energy Tensor. Given our assumptions, the gravita- 
tional 

mm, 

kT * (grav) = ( _ z a_ a + 

2 



tional stress-energy tensor of macroscopic gravity ( - grav ' ) due to connection correlations 
(12) kt; (srav) = {-Z\ vP + ^Q^O"", 



can now be determined, where 

Z a \ivf3 — 2^ a Ate <E i y j a, and Z e ^ = Q^ u . 

The nontrivial components of nTg ( - grav ' ) are 

^tW (s rav ) _ ^( 2 ) fe rav ) _ n 

ki (1) — hi (2) — u 

KT g, (s ™, = 

T (4) (grav) _ _hl 

K1 W ~ %2 



The effective energy density of the gravitational stress-energy tensor of macroscopic gravity 
is 

(grav) rriOi (grav) R 4fl\ 

Kp yti ' = Kip U a U^ = — — > U. 

We note that the sign of p( grav ) is determined by equation (flTT) . The effective energy density 
exterior and interior to a static spherically symmetric matter source have the same functional 
form; in both cases the effective energy density behaves as r~ 2 . The effective pressure of the 
gravitational stress-energy tensor of macroscopic gravity has only a radial component, 

(grav) AH\ (grav) ^ n 

K P radial = —~T = ~ K P^ ' < 

which the macroscopic gravity equations determines to be equal in magnitude, but opposite 
in sign to the effective energy density. If one interprets these results using a fluid formulation, 
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one finds that the gravitational stress-energy tensor of macroscopic gravity can be effectively 
modeled as an anisotropic fluid, with non-trivial radial and zero tangential pressures. Be- 
fore proceeding further, it is worth mentioning that the stress-energy tensor of macroscopic 
gravity satisfies the Weak, Dominant, and Strong energy conditions [15]. 

Coley and Pelavas in a series of two papers [11] have observed something similar; however, 
the process and perspective they followed and what is used in this paper are different. They 
made assumptions about the form of the microscopic geometry and after making a reasonable 
assumption on the nature of the inhomogeneities, explicitly calculated (using Zalaletdinov's 
approach to averaging) the spatially averaged value of the Einstein Field Equations to ex- 
plicitly determine the form of the gravitational correlation tensor. In particular, they started 
with the assumption that the microscopic geometry is described by a spherically symmetric 
metric and determined that the gravitational correlation term behaves as an imperfect fluid. 
Here in this paper, a different perspective is used, a macroscopic perspective; no assumption 
is made about the form of the microscopic geometry, nor any explicit assumptions on the 
nature of the inhomogeneities. The only assumption on the geometry is made at the macro- 
scopic level, and the effects of the inhomogeneities have been incorporated directly into the 
connection correlation tensor. In this paper, the microscopic geometry may be "close" to 
being spherically symmetric and static, in that it will have a group of "almost Killing" vec- 
tors, and only after averaging do these "almost Killing" vectors become true Killing vectors 
of the macroscopic geometry. 

2.5. The Averaged Field Equations. The averaged Einstein Field Equations [21 E] are 

(13) E a /3 = G ae M e p - i^G^M^ = - K (tJ (micro) ) - K T; (grav) 
together with the conservation equation for the averaged microscopic matter 

(14) (tf mici \ a = 0. 

Without assuming anything about the form of the averaged energy density, we obtain the 
following system of differential equations, 

(15) E\ = E\ = e- 2X (----„" + AV - (v') 2 ) = -^ (mat) , 

\r r J 

(16) E% = € -*(^-\-¥)=a!±-kp<"*\ 

l ry /j ry£ ry J qy£i 

(17) E\ = e~^( e ^-\ + 2 -^-)=4^ + Kp^\ 



qy'ii ry 2 ry J ry 2 



and the conservation equation 

(18) p( mat )' = _(p( mat ) _|_p(mat)w_ 
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3. Exterior Solution 

For r > R we have (t^ (micro) ) = 0, (i.e., p (mat) = p (mat) = 0). Equation flUD is identically 
satisfied and equation (TIT]) is easily solved to obtain 

(19) e- 2A W = (1 - 4ht) (l - 2M 

where M is an arbitrary constant of integration. We note that for r > 2M we have the 
constraint hi < 1/4. Subtracting equation ffT7|) from equation ffTB]) we find that 

(20) u' = -A' 
which implies that 

(21) e 2 ^ = Ce~ 2X ^ 

where C is a second arbitrary constant of integration. However, we are able to set C 
to any value through a re-scaling of the coordinate t. Here, for simplicity, we shall set 
C — (1 — 4/ii) -1 and write the macroscopic metric of the exterior of a static spherically 
symmetric source as 

(22) ds 2 = r 2 d9 2 + r 2 sin 2 9 d<f + (1 _^ ^ ( 1 - — ) dr 2 -(l-—)dt 2 

which has the usual coordinate singularity at r = 2M and reduces to the standard Schwarzschild 
solution when hi = 0. Recall, for the standard Schwarzschild metric, the Ricci scalar is 
trivial, however, for the averaged macroscopic geometry determined by the averaged field 
equations and explicitly written as equation fl22l) . we find that the Ricci scalar R = ^ clearly 
showing that the geometry determined by equation f[2"2"j) is different than Schwarzschild. 
It is instructive to rewrite the line element (I2"2l in isotropic coordinates 

x = f 1/vi ~ ihl sin 9 cos 0, 

y = f 1/V1-4/11 Q fa 



z = f 1/vl ~ 4hl cos 9, 
where 




Technically, there is still a two parameter family of connection correlations, but these connection corre- 
lations have no effect on the averaged geometry. 
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which yields a line element of the form 



M 




(23) ds 2 = 1 + i f 2- 2/v /ra7 ( dx 2 + dy 2 + ^ _ 

The macroscopic spacetime described by line-elements (1221) or (|23|) do not approach the 
Minkowski spacetime in the limit as r — > oo or f — > oo (i.e, far field limit), nor should we 
expect it to, since the spacetime is not vacuum, but contains non-trivial corrections from 
averaging. The usual weak-field interpretation is not appropriate here because the weak-field 
limit assumes that the gravitational field is described by Newtonian Gravity, which we have 
assumed is only appropriate on much, much smaller scales. The General Relativistic limit, 
and the consequent weak- field interpretation at infinity only takes place when hi = 0. 



4. Interior Solution 
We begin by defining the mass function m(r) as 



pT 

(24) {l-Ah x )m(r) = - / p (mat) s 2 ds 

2 Jn 



« f r ^ , , \ for < r < R 



R 3 



for R < r 



For r > R we observe that the mass function is a constant, 

tsch 



m(r) 



1 Kp R 3 M s 



[1-Aht) 6 (1-4/n) 



where Ms c h is the usual Schwarzschild mass of the matter source in the absence of connection 
correlations. 

Using this mass function for r < R, equation (|T7|) can be integrated to yield 

^4-<Mli-^ = d-ftl fi-i 



where = 3(1 ~ p 4 Q fel) . Equation (HBJ yields 



(25) p (mat ^ = cie-^ rJ - po 

and where Cx is an arbitrary constant that will be determined by the boundary conditions. 
Substituting equation (|25|) into equation (Tl6|) one obtains a differential equation for v' . 
Integrating, we find 
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where c 2 is another constant of integration. Substituting equation ( 126|) into equation ([25 
we find an expression for the pressure 

— 2 ci + c 2 p (1-^ 



(27) p( mat )(r)=p - 



/ 2 \ V2 

5 ci - c 2Po (l - ^ 



We shall match the interior solution to the exterior solution by requiring that the p( mat \R) = 
0, and that the metric functions, A(r) and z/(r), are continuous at the boundary r = R. This 
is accomplished by setting 



1 A I, ^ 

c 2 = - and ci = p | 1 - ~^ 



1/2 



Therefore, the interior solution of the averaged Einstein Field Equations having constant 
interior energy density po is 



(28) ds 2 = r 2 d6 2 + r 2 sin 2 9 




The continuity condition on the metric functions at the boundary also allows us to interpret 
the the constant of integration M present in the exterior solution. We obtain the condition 

M = = M Sch + -^-M Sch 

1 - 4/ii 1 - 4/i a 6c/l 

which can be now interpreted as the total gravitational mass due to matter together with a 
small correction due to averaging. 

5. Discussion 

The flattening of rotation curves of spiral galaxies is the primary argument for the existence 
of Dark Matter. This Dark Matter is very often modeled as a pressureless dust in a spherical 
halo surrounding the galaxy. However, a pressureless dust yields a density profile that 
behaves as r -3 when in reality, all that is required is a density profile of r~ 2 . It is even more 
curious to note, that the rotational velocity of objects far from the galactic center continues 
to be constant, even though it is commonly assumed that this galactic halo of dark matter has 
some finite boundary, which has yet to be observed. In the context of Macroscopic Gravity, 
we have found that effective energy density due to connection correlations is, p grav oc r~ 2 , and 
that it has the same form both inside and outside the mass generating the gravitational field. 
Consequently, the energy density due to connection correlations can explain the flattening 
of rotation curves of spiral galaxies, without the need to add Dark Matter. A back of 
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the envelope calculation using Newtonian Mechanics^, using a typical observed rotational 
velocity of 200 km/s, for large r, yields hi ~ 10~ 7 . 

It is curious to note that the form of the gravitational stress-energy tensor due to con- 
nection correlations in the simple model presented here takes on the effective form of an 
anisotropic fluid which can explain the flat rotation curves. The possibility of modeling 
dark matter as an anisotropic fluid in spherically symmetric spacetimes has already been 
investigated and more importantly has been shown to be an extremely viable alternative to 
the usual dust models of dark matter [16J. What is new here however, is that Macroscopic 
Gravity now provides the explanation. Clearly, dark matter in the traditional sense is no 
longer required; there is no need to use novel and exotic particles to explain the flat rotation 
curves. Dark matter is the result of averaging general relativity over an appropriate length 
scale, and therefore one should perhaps refrain from using the term "Dark Matter" since it 
may promote confusion. 

There are still many questions that need to be answered in the context of Zalaletdinov's 
theory of Macroscopic Gravity. For example, in order to properly compare the theoretical 
results with astrophysical observations (including the galactic rotation curves mentioned 
here), one must properly study the geodesic equations for both massive and massless test 
particles. A straightforward calculation shows that the geodesic equation for the exterior 
metric given by equation (1221) becomes 

/ \ d 2 u kMc r h 9 

(29) — + u = —jp: + Mdsch^ + 4M 

Newton GR correction MG correction 

where u = 1/r, Mgch is the Schwarzschild mass, L is the angular momentum, k = for 
massless particles, and k = 1 for massive particles. The form of the geodesic equation is 
nearly the same as in the traditional version of the geodesic equation for the Schwarzschild 
metric. We note the appearance of an additional term linear in u that is not present in 
the geodesic equation for the standard Schwarzschild spacetime. Through a re-scaling of 
the independent variable (0 = \Jl — Ahi(f>), one can show that the geodesic equation ( l29l) 
becomes 

d 2 u kM ~ 2 

+ u = — - + 3Mu 2 
d<P 2 L 2 

Since M ~ Ms c h, we find that there would only be a small change in the astrophysical ob- 
servations. However, any change to the geodesic equations ought to be thoroughly analyzed. 

The analysis of gravitational and cosmological models based on Zalaletdinov's Macroscopic 
gravity is by no means complete. Even in the spherically symmetric and static case studied 
here, the effects of the electric part of the connection correlation tensor Z is yet unknown. 



2 For now this is just a toy model, the assumption that Newtonian mechanics is adequate for the description 
of the dynamics of galaxies may be questionable, as the equations for geodesies also change because of the 
averaging. 
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The electric part of the connection correlation tensor may significantly alter or on the other 
hand may fine tune the solution to the galactic rotation curve problem. The electric part 
of the connection correlation tensor may also provide a natural explanation for the accel- 
eration of the universe in a cosmological context. Furthermore, the assumption that the 
connection correlation tensor Z is invariant under the same Lie group of Killing vectors as 
the macroscopic metric, although ensuring compatibility between them, is restrictive. It will 
be of great interest to relax this invariance condition and study the effects of more general 
connection correlations even in the static spherically symmetric case studied here. 

Nevertheless, the result presented in this paper illustrates that Zalaletdinov's Macroscopic 
Gravity provides an extremely promising explanation for the flattening of rotation curves 
and therefore this averaged theory of gravity merits further investigation. 

6. Acknowledgements 

The author would like to thank Roustam Zalaletdinov for his tremendous patience, under- 
standing, and guidance. The author would also like to thank Alan Coley for his comments. 
MG Eqs Code, [T7|, GRTensorII[IS] and the computer algebra package Maple were used 
extensively for the calculations present in this paper. This work was supported through a 
grant from the Natural Sciences and Engineering Research Council (NSERC) of Canada. 

References 

[1] C. M. Will, Theory and Experiment in Gravitational Physics (Revised Edition), (Cambridge University 
Press, Cambridge, 1993). 

[2] S. Perlmutter et al., Astrophysical Journal, 517, 565 (1999); A. G. Reiss et al., Astrophysical Journal, 
116, 1009, (1998). 

[3] V. C. Rubin, Science, 220, 1339, (1983); V. C. Rubin, W. K. Ford, and N. Thonnard, Astrophysical 

Journal, 238, 471, (1980); V. C. Rubin, W. K. Ford, N. Thonnard, and D. Burstein, Astrophysical 

Journal, 261, 439, (1982); V.C Rubin, D. Burstein, W. K. Ford, and N. Thonnard, Astrophysical 

Journal, 289, 81, (1985); Y. Sofue and V. C. Rubin, Ann. Revs. Astron. and Astrophys., 39, 137, (2001); 

A. Borriello and P. Salucci, Monthly Notices of the Royal Astronomical Society, 323, 285 (2001); M. 

Persic, P. Salucci and F. Stel, Monthly Notices of the Royal Astronomical Society, 281, 27, (1996); L. 

Bergstrom, Rep. Prog. Phys., 63, 793 (2000); J. P. Ostriker and P.J.E. Peebles, Astrophysical Journal, 

186, 467, (1973); S. Tremaine, Physics Today, 45, 28-36, (1992). 
[4] G. F. R. Ellis, in General Relativity and Gravitation, ed. B. Bcrtotti, F. de Felice and A. Pascolini 

(Dordrecht: Reidel, 1984), 215; G. F. R. Ellis, W. Stoeger, Class. Quantum Grav., 4, 1697, (1987). 
[5] R. M. Zalaletdinov, Gen. Ret Grav., 24, 1015-1031, (1992); R. M. Zalaletdinov, Gen. Rel. Grav., 25, 

673, (1993); R. M. Zalaletdinov, Bull. Astron. Soc. India, 25, 401, (1997). 
[6] R. M. Zalaletdinov, Int. Jour. Mod. Phys. A, to appear (2008), preprint arXiv:gr-qc/0801.3256; 

R. M. Zalaletdinov, "Averaging Problem in Cosmology and Macroscopic Gravity" , preprint 

|arXiv:gr-qc/0701116| 

[7] A. Krasihski, Inhomogeneous Cosmological Models, (Cambridge University Press, Cambridge, 1997) 
[8] M. F. Shirokov, I. Z. Fisher, Sov. Astron. J., 6, 699 (1963). Reprinted in: Gen. Rel. Grav. 30 1411- 
1427, (1998); T. W. Noonan, Gen. Rel. Grav., 16, 1103, (1984); T. W. Noonan, Gen. Rel. Grav., 17, 
535, (1985); R. A. Isaacson, Phys. Rev., 166, 1263, (1968); R. A. Isaacson, Phys. Rev., 166, 1272, 
(1968); T. Futamase, Phys. Rev. Lett, 61, 2175, (1988); T. Futamase, Phys. Rev. D, 53, 681, (1996); 



14 



R. J. VAN DEN HOOGEN 



S. Bildhauer and T. Futamase, Gen. Rel. Grav., 23, 1251-1264, (1991); N. V. Zotov and W. Stoeger, 
Class. Quantum Grav., 9, 1023 (1992); M. Kasai, Phys. Rev. D, 47, 3214, (1993); J. P. Boersma, Phys. 
Rev. D, 57 (1998), 798; T. Buchert, Gen. Rel. Grav., 32, 105-125 (2000); T. Buchert, Gen. Rel. Grav., 
33, 1381-1405, (2001). 

[9] A. Paranjape and T.P. Singh, Phys. Rev. D 76, 044006 (2007); A. Paranjape, "A Covaria nt Road to Spa- 
tial Averaging in Cosmology: Scalar Corrections to the Cosmological Equations" Preprint gr-qc/0705238 
(2007); 

[10] A. A. Coley, N. Pelavas, and R. M. Zalaletdinov, Phys. Rev. Lett, 95, 151102, (2005); 
[11] A. A. Coley and N. Pelavas, Phys. Rev. D, 74, 087301, (2006); A. A. Coley and N. Pelavas, Phys. Rev. 
D, 75, 043506, (2007). 

[12] J. Plebahski and A. Krasihki, An Introduction to General Relativity and Cosmology, (Cambridge Uni- 
versity Press, Cambridge, 2006). 
[13] R. J. van den Hoogen and R. M. Zalaletdinov, (2006) [in preparation]; 
[14] R. M. Zalaletdinov, private communication. 

[15] S.W. Hawking and G.F.R. Ellis, The Large Scale Structure of Space-time, (Cambridge University Press, 
Cambridge, 1973). 

[16] K. Lake Phys. Rev. Lett., 92, 051101, (2004); S. Bharadwaj and S Kar, Phys. Rev. D, 68, 023516, 
(2003); U. Nucamendi, M. Salgado, and D. Sudarsky, Phys. Rev. D, 63, 125016, (2001). 

[17] R. J. van den Hoogen and R. Zalaletdinov, The MG Eqs Code [A computer algebra code for generation 
and manipulation with the equations of Macroscopic Gravity] (St. Francis Xavier University, 2007). 



[18] GRTensorll is a software package that runs within MAPLE. It is distributed freely at http: / / grtensor.org 



R. J. van den Hoogen, Department of Mathematics, Statistics, and Computer Science, St. 
Francis Xavier University,, Antigonish, Nova Scotia, Canada, B2G 2W5 
E-mail address: rvandenh@stfx.ca 



